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ABSTRACT 

This  report  documents,  for  the  first  time,  the  wavenumber  and  frequency 
characteristics  of  the  mean-shear  (MS)  and  turbulence-turbulence  (TT)  par¬ 
tial  wall  pressures.  Results  show  that  both  the  TT  and  MS  pressures  are 
important  for  the  convective  pressures,  but  that  the  TT  pressure  dominates 
in  the  subconvective  wavenumber  range.  At  low  frequencies  both  the  buffer 
layer  and  logarithmic  region  contribute  significantly  to  the  MS  convective 
pressure;  at  higher  frequencies  the  buffer  layer  is  dominant.  The  partial 
pressures  were  obtained  from  two  large  eddy  simulations  (LES)  of  a  fully- 
developed  turbulent  channel  flow;  they  were  computed  from  the  TT  and 
MS  source  terms  from  four  regions  of  the  channel  (viscous  shear  layer, buffer 
layer,  logarithmic  region  and  upper  channel).  The  Reynolds  number  for  the 
simulations  was,  Reg*  =  Ur5*lv  pa  24  (Rcj  =  171  j;  the  non-linear  terms 
in  the  simulation,  as  well  as  those  pressure  source  terms  were  dealiased  us¬ 
ing  the  2/3  rule.  The  first  simulation,  LESS,  used  the  rotational  form  of 
the  nonlinear  terms,  whereas  the  second  simulation,  LES9,  used  the  skew- 
symmetric  form.  Significant  differences  between  the  two  simulations  appear 
only  in  the  subconvective  wavenumber  range  for  the  TT  spectra  that  has 
been  integrated  over  the  entire  channel,  pointing  to  the  fact  that  the  skew- 
symmetric  form,  while  more  expensive,  computationally,  is  more  accurate. 

EXECUTIVE  SUMMARY 

This  work  is  part  of  an  effort  to  develop  a  hybrid  RANS/ Analytical/ Statistical  nu¬ 
merical  tool  for  the  prediction  of  low  Mach  number,  hydrodynamically-generated  soimd 
in  the  absence  of  cavitation.  This  tool  is  based  on  a  mathematical  formulation  for  the 
mean-shear  (MS)  component  of  the  pressure  in  which  RANS  simulations  provide  the 
mean  flow  variables,  while  models  based  on  turbulence  simulations  provide  the  neces¬ 
sary  information  regarding  the  fluctuating  somrces  and  their  wall-normal  correlations. 
Using  turbulence  simulations,  all  the  elements  of  modeling  can  be  tested  since  the  ex¬ 
act  three-dimensional  source  terms  and  their  resultant  wall  pressmes  can  be  obtained. 
In  the  work  reported  herein  we  characterize  the  wall-pressure  fluctuations  due  to  the 
MS  and  TT  pressures  in  both  the  convective  and  subconvective  wavenumber-frequency 
ranges.  This  work  is  novel  in  three  respects:  (1)  it  characterizes  the  wavenumber- 
frequency  characteristics  of  the  individual  MS  and  TT  pressures,  (2)  it  characterizes 
these  pressures  in  the  subconvective  wavenumber-frequency  range  and  (3)  the  MS  and 
TT  pressures  have  been  decomposed  into  their  components  from  four  layers  parallel  to 
the  wall.  Thus,  this  work  provides  a  first  glimpse  into  the  nature  of  the  MS  and  TT 
pressures  in  a  form  appropriate  for  hydroacoustic  modeling  and  it  provides  a  database 
which  can  be  used  for  validating  models  of  the  MS  (and  TT)  pressures;  in  addition 
it  shows  the  relationships  between  the  the  MS  and  TT  velocity-field  sources  and  their 
pressures  in  wavenumber-frequency  space  —  this  will  help  with  physics-based  mod- 
ehng.  The  wall-pressure  and  velocity  field  databases  were  obtained  from  LES  of  a 
fully-developed,  tiubulent  channel  flow  with  Reg*  =  UrS*lu  «  24  (Re^  =  171). 
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INTRODUCTION 

The  US  Navy  is  interested  in  tm-bulent  wall  pressure  fluctuations  (WPF)  because 
they  are  an  important  source  of  energy  for  flow-induced  noise.  The  low- wavenumber, 
high-frequency  (with  respect  to  the  convective  ridge)  components  (or  subconvective) 
of  WPFs  provide  the  forcing  function  for  the  structures  underlying  the  flow.  The  con¬ 
vective  components  provide  the  forcing  function  for  flow-induced  noise  from  trailing 
edge  flows.  The  nature  of  the  subconvective  pressure  has  been  difficult  to  determine, 
due  to  its  low  energy  (with  respect  to  the  convective  ridge),  whereas  the  nature  of  the 
pressure  at  the  convective  ridge  has  been  easier  to  determine  due  to  its  higher  energy 
levels.  Nonetheless,  for  modelling  purposes,  enough  information  from  experiments  has 
been  hard  to  come  by  due  to  lack  of  three-dimensional  data  with  which  to  characterize 
the  velocity  field  wall  pressure  somrces.  Chang,  Piomelli  and  Blake^,  using  a  velocity 
field  database  from  a  low- Re  turbulent  channel  flow  (TCF)  decomposed  the  WPFs  into 
their  mean-shear  (MS)  and  turbulence-turbulence  (TT)  components  (partial  pressures) 
and  further,  due  to  four  regions  of  the  channel  (viscous  shear  layer,  buffer  layer,  log¬ 
arithmic  region  and  upper  channel).  They  showed  in  one-dimensional  streamwise  and 
spanwise  wavenumber  space  which  terms  and  regions  are  responsible  for  the  dominant 
partial  pressures.  However,  the  one-dimensional  spectra  show  essentially  the  character 
of  the  peak  of  the  convective  ridge,  disguising  the  nature  of  the  partial  pressures  away 
from  the  convective  ridge,  particularly  in  the  subconvective  wavenumber  range.  Chang, 
Piomelli  and  Blake^  showed  that  the  spectra  of  the  individual  TT  pressures,  ny(Rr) 
(where  Rr  denotes  one  of  the  foiu  regions  of  the  channel)  have  two  or  more  order  of 
magnitude  greater  energy  in  the  lowest  wavenumber  modes  than  the  TT  pressure  due 
to  all  terms  and  regions.  There  is  speculation  that  this  high  energy  is  due  to  low- 
wavenumber  cross-channel  coherence.  Studying  the  kx—oj  characteristics  of  the  partial 
pressures  may  shed  some  more  light  on  the  natme  of  the  low-wavenumber  behavior. 

In  Chang^  modelling  the  convective  ridge  pressure  was  addressed  through  the  eval¬ 
uation  of  a  model  for  the  MS  pressure.  A  complete  kx—u  spectrum  of  the  MS  pressure 
is  necessary  to  guide  model  development,  provide  data  for  statistical  functions  and 
provide  data  for  validation. 

The  objective  of  this  work  was  to  study  the  wavenumber-frequency  characteristics 
of  the  partial  pressmes,  particularly  to  assess  the  subconvective  behavior  and  to  study 
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computational  issues  that  affect  it.  Two  long  time-series  of  partial  pressures  was  ob¬ 
tained  from  a  LES  of  a  TCP  simulation.  FVom  this  database,  the  wavemunber-frequency 
spectra  of  the  partial  pressures  was  computed  and  evaluated. 


MATHEMATICAL  FORMULATION 

In  this  section  we  show  the  various  forms  for  the  nonlinear  terms  and  suggest  how 
they  might  impact  computation  of  the  wall  pressure.  The  convection  form,  as  derived 
from  integral  calculus,  is  given  by 

h  =  uVu,  (1) 


the  divergence  form  of  h  is  given  by 


h  =  V  •  (uu) , 


(2) 


the  rotational  form  is  given  by 

h  =  u  X  a;  -t-  ^  V  (u  •  u)  (3) 

and  the  skew-symmetric  form  is  given  by 


h=i[u-(Vu)-HV-(uu)].  (4) 

The  form  of  the  nonlinear  terms,  h,  in  the  Navier-Stokes  equations  may  be  impor¬ 
tant  for  the  computation  of  the  wall  pressure  directly  and  indirectly.  The  form  of  the 
nonlinear  terms  in  the  Poisson  equation  may  have  a  direct  affect  on  the  wall  pressure. 
The  TT  source  terms  are  nonlinear  velocity  products  derived  directly  from  h.  The 
product  of  two  sinusoidal  signals  are  themselves  signals  that  have  wavenumbers  which 
are  the  sum  and  difference  of  the  two  signal’s  respective  wavenumbers,  e.g.. 


sin  a  cos  /3 


^  cos(q:  -  ^  cos(a  -f  /?). 


(5) 


The  nonlinear  interactions  result  in  the  transfer  of  energy  primarily  into  modes  at  the 
convective  ridge.  Energy  is  transferred  to  other  parts  of  the  spectrum,  including  the 
subconvective  wavenumber  range.  The  subconvective  wavemunber  range  has  energy 
levels  that  are  at  least  20  dB  down  from  the  convective  ridge  so  it  is  conceivable  that 
sTTiflll  changes  in  how  h  are  computed  can  affect  the  subconvective  wavenumber  range 
characteristics.  For  simplicity  the  divergence  form  has  always  been  used  for  computing 
the  nonlinear  source  terms. 

The  pressure  is  affected  indirectly,  in  the  computation  of  the  velocity  fields  them¬ 
selves  in  the  N-S  equations.  When  discretized,  the  convection  form  does  not  conserve 
linear  momentum  or  kinetic  energy;  the  divergence  form  conserves  only  linear  momen¬ 
tum.  This  leads  to  numerical  instabilities  in  tmbulence  simulations.  However,  it  has 
been  shown  (see  Zang^)  that  switching  between  the  convective  and  divergence  forms 
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(the  “alternating”  form)  results  in  a  well-behaved  simulation.  On  the  other  hand,  the 
rotation  and  skew-symmetric  forms  conserve  both  linear  momentum  and  kinetic  en¬ 
ergy.  It  can  be  seen  that  the  skew-symmetric  form  is  the  average  of  the  divergence 
and  convection  forms,  and  thus  closely  related  to  the  alternating  form.  The  rotation 
form  only  requires  the  computation  of  six  derivatives  while  the  skew-symmetric  form 
requires  18;  the  alternating  form  would  take  six  and  12  derivative  computations  on 
alternating  time  steps.  For  efficiency  reasons,  the  rotation  form  is  the  most  popu¬ 
lar,  used  for  instance  in  Kim  et  al.,^  and  Choi  and  Moin^.  Zang^  showed  that  for  a 
fully-developed  turbulent  channel  flow  the  form  of  the  nonlinear  terms  made  differ¬ 
ences  in  the  high-wavenumber  tail  of  the  velocity  spectra  but  did  not  affect  the  lower 
wavenumbers.  However,  in  a  simulation  of  homogenous  turbulence  it  was  shown  that 
the  rotation  form  significantly  underpredicts  the  turbulence  kinetic  energy  and  velocity 
spectra.  Also  compared  were  the  integral  scales,  microscales,  skewness  and  flatness;  it 
was  found  that  the  skew-symmetric  and  alternating  forms  had  results  much  closer  to 
the  Galerkin  method*  results  than  did  the  rotational  results. 

When  the  rotational  form  is  used  the  resolved  kinetic  energy  term  is  incorporated 
into  the  pressure  term  such  that  the  computed  pressme  is  P  =  p  -f  |u  •  u,  where  p  is 
the  dynamic  pressure.  Furthermore,  the  dynamic  subgrid-scale  model  only  computes 
the  anisotropic  part  of  the  SGS  stresses  and  the  isotropic  part  is  assumed  to  be  added 
to  the  pressure  term.  For  the  present  work  the  wall  pressure  is  computed  directly  from 
the  velocity  fields.  The  nonlinear  terms  were  evaluated  with  the  rotational  form  for 
LESS  and  the  skew-symmetric  form  for  LES9. 

In  both  cases,  the  nonlinear  terms  were  dealiased  in  the  x  and  2  directions  using 
the  2/3  rule.  In  this  method,  the  time  advancement  is  performed  on  a  “3/2-grid”  in 
the  X  and  z  directions;  the  velocity  fields  are  then  converted  to  Fourier  space  where 
the  top  1/3  modes  in  each  direction  are  truncated,  removing  the  aliasing  errors  from 
the  velocities;  the  velocity  fields  are  then  transformed  back  into  physical  space.  Restart 
files  were  written  out  on  the  3/2  grid;  these  files  were  Fourier  transformed,  the  top  1/3 
modes  were  truncated  (they  are  zero,  anyway),  inverse  Fourier  transformed  and  then 
written  to  file.  For  reading  into  the  simulation,  the  velocity  restart  files  were  spectrally 
expanded  onto  the  3/2  grid. 

Computation  of  the  frequency  spectra  entails  obtaining  long  time  series  of  pressure 
data.  The  LES  was  run  on  a  SGI  Origin  200  (RIOOOO  at  188  MHz).  The  CFL  number 
was  set  to  0.14324  which  corresponds  to  a  time  step  AtUe/5  =  0.025.  This  resulted 
in  maximum  Courant  numbers  of  approximately  0.30,  well  below  the  stability  limit  of 
0.50. 

The  LES8  database  consists  of  51  files  with  64  realizations  per  file  for  a  total  of  3264 
realizations  spanning  a  time  of  3265 /Ug  {9.865 /ut)-  For  each  realization,  10  partial 
pressures  were  saved:  four  fields  each  of  TT  and  MS  pressure  corresponding  to  the  four 
regions  and  the  TT  and  MS  pressures  from  all  regions. 

*In  the  Galerkin  method  the  nonlinear  terms  are  computed  in  wavespace  rather  than  in  physical 
space.  This  method  is  considered  more  accurate,  but  is  much  more  computationally  intensive. 
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The  LES9  database  consists  of  20  files  with  64  realizations  per  file  for  a  total  of  1280 
realizations  spanning  a  time  of  V2&5/Ut  (9.86d/«T).  For  each  realization,  11  partial 
pressures  were  saved:  fom:  fields  each  of  TT  and  MS  pressure  corresponding  to  the 
four  regions,  the  TT  and  MS  pressures  from  all  regions  and  the  total  pressure  from  all 
regions. 

Computation  of  partial  pressures 


Table  1.  Regions  of  the  chaimel. 


Region 

Limits 

Description 

1 

0  <  <  5 

Viscous  shear-layer 

2 

5  <  y+  <  30 

Buffer  layer 

3 

30  <  y+  <  180 

Logarithmic  region 

4 

180  <  y+  <  360 

Upper  channel 

Since  the  wall-pressure  depends  on  sources  distributed  throughout  the  entire  flow 
domain  it  is  difficult  to  ascertain  where  in  the  boundary  layer  the  dominant  soiurces 
are  located.  One  advantage  of  turbulence  simulations  is  that  they  generate  three- 
dimensional  velocity  fields  from  which  all  the  somce  terms  can  be  computed.  In  this 
investigation  the  product  of  the  somces  by  the  Green’s  function,  which  is  the  true 
measure  of  the  influence  of  the  sources  on  the  wall-pressiue,  has  been  computed.  The 
channel  was  divided  into  four  regions  (Table  1)  roughly  corresponding  to  the  viscous 
shear  layer,  buffer  and  logarithmic  regions,  and  the  upper  part  of  the  channel.  The 
contribution  of  the  source  terms  was  integrated  over  each  region.  It  is  assiuned  that  the 
partial  pressures  are  evaluated  at  the  lower  channel  wall,  even  though  symmetry  has 
been  used  in  the  data  reduction  and  partial  pressmres  have  been  averaged  over  both  the 
upper  and  lower  channel  walls.  The  lower  wall  will  be  referred  to  as  the  “near-wall,” 
the  upper  wall,  the  “far-wall.” 

The  pressure  source  terms  were  computed  in  physical  space  on  the  3/2  grid. 
Nine  terms  were  computed:  the  six  TT  terms,  their  sum,  the  MS  term  and  the  to¬ 
tal.  They  were  then  Fourier-transformed  to  yield  their  complex  coefficients  in  planar 
wavenumber  space,  Tij(K,  y).  The  integral  of  the  product  of  the  source  terms  and 
Green’s  function 

Pij{K,  Rr)  =  /  7;,(K,  r))G{K,  y  = -1,  r))dr).  (6) 

(where  771  (r)  and  772  (r)  are  the  lower  and  upper  integration  limits  for  each  region,  and  r 
denotes  the  vertical  regions  described  in  Table  1)  was  then  computed  using  Chebychev 
integration.  The  complex  partial  pressme  spectrum  due  to  the  source  term  Tij  from 
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region  r  is  denoted  by  p(K,  Rr).  The  complex  spectra  for  the  various  terms  and  regions 
could  then  be  combined;  e.g.,  the  partial  pressure  due  to  T12  from  regions  2,  3  and  4  is 

g^{K,R234)=p'S'{K.It2)+pS’(K,m)+p'S'(K,R4).  (7) 

The  total  TT  partial  pressure  is  the  summation  of  all  the  TT  partial  pressures, 

Rr)  =  (S) 

i=l  j=l 

and  the  total  partial  pressure  is 

Rr)  =  Rr)  +  p^{K,  Rr).  (9) 

The  total  pressure  is  i?1234). 

The  magnitudes  of  the  complex  spectra  were  computed  by 


n(K,  Rr)  =  plj{K,  Rr)pij{K,  Rr)  (no  summation)  (10) 

(where  f*  is  the  complex  conjugate  of  /),  and  averaged  over  multiple  realizations.  It 
should  be  noted  that  the  partial  pressures  are  defined  by  the  integral  over  all  frequencies, 


roo 

U{K,Rr)=  U{K,u,Rr)dw 

J  —  oo 

(11) 

where  n(K,a;, /?r)  is  the  partial  pressure  in  four-space, 
were  obtained  by  the  integrations 

The  one-dimensional  spectra 

roo 

TT{kx,Rr)=  /  ll(K,Rr)dkz 

(12) 

and 

TOO 

n{kz,Rr)=  /  U{K,Rr)dkx 

J — oo 

(13) 

The  one-dimensional  spectra  were  normalized  by  t^S. 


Computation  of  partial  pressure  spectra 

The  wavevector-frequency  spectra  shown  in  this  Section  were  computed  by  the  fol¬ 
lowing  methodology:  the  data  were  divided  into  Nt  =  384  point  segments  spanning 
the  period,  T,  with  50  percent  overlap.  This  gave  a  reasonable  compromise  between 
number  of  spectra  realizations  and  frequency  resolution  (the  frequency  resolution  is 
approximately  Aur/S  and  the  0.16 /ur  sampling  rate  gives  a  maximum  frequency  of  ap¬ 
proximately  760ur/5).  A  Hanning  window  was  used  on  the  time  segments,  while  no 
windowing  was  necessary  in  the  x  and  ^  directions  in  which  the  flow  is  periodic.  The 
Hanning  window  was  applied  in  time  to  the  total,  resolved  pressrue,  p,  normalized  by 
the  long  time  average  value  of  r^,. 


p{Xi,Zj,ti)  =p{Xi,Zj,ti)  - 


(14) 
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f* 


where  i  =  \...Nx,  j  =  l...Nz,  and  I  =  A  three-dimensional  Foiuier 

wavenumber-frequency  transform  of  each  time  segment  was  taken  and  from  that,  an 
averaged  spectral  density  was  obtained  from  the  mean  of  the  spectral  densities  from 
the  m  overlapping  segments 

1  m 

$  (iAk,,  }AK,  lAu)  =  -  E  (P  P*).  •  (15) 

^  S=1 

Here  p  denote  the  Fourier  coefficients  of  the  pressure  and  p*  are  their  complex  conju¬ 
gates.  The  three-dimensional  spectral  density  should  satisfy  Parseval’s  theorem  which 
states  that  the  “mean  square  value  of  a  periodic  function  is  equal  to  the  sum  of  squares 
of  all  it  Fourier  coefficients,®” 


-1  Nt  Nz  Nz  Nt  Nz  Nx 

WWIV  S  E  =  E  E  E  ^ 


(16) 


where  the  left  hand  side  is  the  mean-square  (MS)  pressure.  Since  the  pressure  has  been 
windowed  in  time,  its  spectral  density  does  not  satisfy  (16)  exactly;  is  then  normalized 
to  satisfy  (16).  The  two-dimensional  kx—u)  spectra  is  given  by  the  summation  over  all 
the  spanwise  wavenumber  modes, 


Nz 

^iAkxJAuj)  =  'Z-^.  (17) 

J=1 

The  two-dimensional  spectral  density  extends  over  I  =  1.  ..Nt  frequency  bins  and 
i  =  l...Nx  wavenumber  bins  and  has  the  symmetry  property 

^iAkx,  lAu)  =  ^[{Nx-i  +  2)Akx,  {Nt-l  +  2)  Aw]  (18) 

In  the  kx  —  oJ  plots  the  wavenumber  bins  Nx/2+2  <i  <  Nx  have  been  translated  to  the 
bins  —Nx/2  +  l<i<—l;  with  this  convention,  a  wave  that  is  convecting  downstream 
will  have  either  w  or  kx  negative.  In  the  two-dimensional  spectra  plots,  the  convective 
ridge,  representing  energy  convecting  downstream,  is  shown  in  the  negative  kx,  positive 
u  quadrant.  Constant  frequency  cuts  of  streamwise  wavenumber  spectra,  extracted 
from  this  two-dimensional  spectra,  show  just  the  negative  kx,  i.e.,  where  the  convective 
ridge  lies. 

The  one-dimensional  point  frequency  spectra  is  obtained  from  the  double  summation 

<l>{lAu)  =  E  E  ^  {i^kx,jAkx,  lAuj) ,  (19) 

j=l  i=l 

and  the  one-dimensional  wavenumber  spectra  are  given  by 
Nz  Nt 

(j){iAkx)  =  E  E  ^  {iAkx,jAk;„  lAu) ,  (20) 

j=i  1=1 
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and 


^O'AS:,)  =  E  E  *  ( Aa.) .  (21) 

i=l  J=1 

These  are  functions  which  are  symmetric  about  their  folding  frequency  (or  wavenum¬ 
ber);  they  have  been  converted  to  one-sided  spectra  by  multiplying  all  components, 
except  the  zero  frequency  (or  wavenumber)  component,  by  a  factor  of  2,  thus  preserv¬ 
ing  Parseval’s  identity. 

Note  that  two  normalizations  have  been  used  in  our  previous  works: 

1.  The  wall-pressure  spectra  shown  in  Chang  et  al7  and  Chang^  are  derived  from 
three-dimensional  kx  —  kz  —  u  spectra.  In  these  works,  three-dimensional  spectra 
was  normalized  using  a  Parseval’s  equation  in  which  the  integral  of  the  spectra 
equalled  the  m.s.  of  the  input  signal.  The  various  summations  e.g.,  (19)-(21) 
were  normalized  such  that  their  integrals  are  equal  to  the  m.s..  The  rationale 
for  this  approach  was  for  dimensional  consistency:  it  is  assumed  that  the  three- 
dimensional  spectra  has  the  units  of  p^L'^T ;  integrating  over  kz  (summing  all  kz 
components  and  multiplying  by  Akz)  gives  the  units  p^LT;  subsequent  integration 
over  kx  to  obtain  the  point  frequency  spectrum  gives  the  units  p^T  which  is 
typically  found  in  the  literature  (see  e.g.,  Keith  et  al.^). 

2.  The  second  normalization,  used  in  the  present  work,  as  well  as  for  the  partial 
pressure  spectra  in  Chang^  is  based  on  the  principle  that  the  summation  of  the 
spectra  equals  the  m.s..  In  this  case  there  is  no  explicit  way  to  change  the  units 
when  performing  the  various  integrations.  It  will  be  assumed  that  the  units  of  the 
kx  —  kz  —  uj  spectra  are p^L^T ;  the  kx  —  kz  spectra p^L^\  kx—io  and  kz  —  uj  spectra 
p^LT ;  one-dimensional  kx  and  kz  spectra  p^L]  and  one-dimensional  u  spectra  p^T. 

To  convert  from  the  first  to  second  normalization,  simply  multiply  the  streamwise 
wavenumber  spectra  by  Akx,  the  spanwise  wavenumber  spectra  by  Akx  and  the  fre¬ 
quency  spectra  by  Au. 


PRESENTATION  OF  RESULTS 


Mean  velocity  and  turbulence  intensity  profiles 

In  this  subsection  we  show  the  statistics  of  the  partial  pressures.  The  mean  ve¬ 
locity  profiles  are  shown  in  Figure  1.  Comparison  to  DNS  results  show  that  they  are 
well-resolved  and  have  the  expected  behavior  in  the  various  regions  of  the  channel. 
Turbulence  intensity  profiles  are  shown  in  Figure  2. 

Mean  square  variation  of  pressure  with  time 

Figure  3  shows  the  variation  of  the  TT,  MS  and  total  pressures  due  to  sources 
over  the  entire  channel[R(1234)].  The  m.s.  pressures  have  been  normalized  by  (pu^S)^, 
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Figure  1.  Viscous-scaled  mean  velocity  profiles;  —  :  LESS; - ;  DNSS. 

while  time  has  been  normalized  by  5/Ui.  It  shows  clearly  that  the  MS  pressure  is  almost 
always  greater  than  the  TT  pressure  and  that  the  largest  variations  in  the  total  pressure 
are  highly  correlated  with  those  variations  in  the  MS  pressure.  By  contrast,  the  TT 
pressure  has  a  much  smaller  variation  than  either  the  MS  or  total  pressure  (even  though 
it  can  be  seen  that  the  TT  variation  are  correlated  with  the  MS  and  total  pressure). 
This  result  is  an  indication  that,  at  least  in  a  broad-band  sense,  and  over  the  entire 
domain,  high  amplitude  events  are  more  closely  associated  with  the  MS  pressure  than 
the  TT  pressure.  This  is  consistent  with  the  findings  of  Johansson  et  al}^. 

Point  frequency  spectra 

The  objective  of  this  subsection  is  to  show  that  the  LESS  compares  very  well  with 
previous  simulations  and  recent  experiments. 

Figure  4  compares  the  point  frequency  spectra  for  LESS,  DNS5,  and  the  TBL  sim¬ 
ulation  of  Singer^^.  Also  shown  are  two  recent  experimental  results  by  Farabee  and 
Casarella^^  and  Gravante  et  with  Re-r  values  similar  to  that  of  Singer.  The  impor¬ 
tant  parameters  for  the  point  frequency  spectra  are  listed  in  Table  2.  Figure  4(a),  the 
spectra  normalized  by  outer  variables,  shows  that  for  the  low  frequencies,  the  LES  and 
DNS  spectra  are  identical.  For  coS/ur  >  100  the  LES  and  DNS  spectra  diverge  rapidly, 
due  to  the  simulations  respective  wavenumber  resolutions  —  the  point  frequency  spec¬ 
tra  is  the  integration  over  all  wavenumbers,  and  it  can  be  seen  in  Figure  4,  that  as 
the  frequency  increases,  the  amoxmt  of  energy  truncated  by  the  wavenumber  cutoflF  in¬ 
creases.  In  fact,  at  about  uj5/ut  =  140,  only  about  one-half  the  convective  ridge  is 
resolved  and  the  difference  between  the  LES  and  DNS  is  about  3  dB. 


9 


Figure  2.  Turbulence  intensity  profiles;  —  :  LESS; -  :  Kim  et  al^; 

o  :  Kreplin  and  Eckelmann^ 


The  simulation  of  Singer  has  lower  energy  than  LESS  and  DNSS  in  the  lower  fre¬ 
quencies,  but  is  virtually  identical  with  that  of  Gravante  et  al}^.  The  data  of  Farabee 
and  Casarella^^  is  slightly  lower.  This  variation  in  the  low  frequencies  is  consistent  with 
the  trends  shown  in  Keith  et  who  compared  point  frequency  spectra  over  a  wide 
range  of  Re.  Farabee  and  Casarella^^,  on  the  other  hand,  showed  that  with  a  two-fold 
Variation  in  Re,  this  frequency  range,  that  includes  the  peak  of  the  spectra,  should 
collapse  to  a  single  curve. 

It  is  apparent  that  the  higher  Re  spectra  has  a  much  larger  spectra  peak  region  than 
does  the  lower  Re  LESS  and  DNS5  data.  This  is  consistent  with  the  notion  that  higher 
Re  flows  have  more  energy  in  the  scales  small  compared  with  S.  The  overlap  range,  as 
proposed  by  Bradshaw^"*  should  have  a  uj~^  decay  rate.  The  spectra  of  Gravante  et  al.^^ 
is  the  only  one  that  exhibits  an  extended  frequency  range  with  such  behavior.  The  TCF 
data  has  no  such  region,  since  low  Re  flows  have  essentially  no  wavenumber  separation 
between  the  larger,  energy  producing  scales  and  the  smaller  energy  dissipating  scales. 
The  data  of  Singer  may  have  a  small  frequency  range  that  has  a  decay  rate;  however, 
it  too,  is  attenuated  rapidly  by  the  simulation’s  wavenumber  resolution.  The  data  of 
Farabee  and  Casarella^^  has  an  approximate  decay  rate,  but  no  u}~^  region. 

The  inner-scaled  spectra.  Figure  4(b),  shows  that  inner  scaling  does  a  fairly  good 
job  of  collapsing  all  of  the  spectra,  with  differences  explainable  by  wavenumber  reso¬ 
lution  arguments  (values  of  u;^,  the  frequencies  where  one-half  of  the  convective  ridge 
is  resolved,  are  listed  in  Table  2).  In  the  range  tt;'*'  <  1  (a;"*"  =  ojS'^/ut)  the  TCF  data 
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Figure  3.  Temporal  variation  of  m.s.  (mean-square)  TT,  MS  and  total  pres¬ 
sures  (R1234); - :  TT, - - :  MS,  and  —  :  total. 


and  Gravante  et  are  identical.  Both  the  DNS  and  Gravante  et  have  a  range 
with  a  decay  rate;  the  DNS  data  diverges  at  about  >  2  due  to  wavenumber 
resolution,  whereas  the  data  of  Gravante  et  al}^,  has  a  behavior  out  to  its  largest 
frequency;  this  is  because  it  has  =  5.  As  expected,  the  spectra  of  Singer  decays 
most  rapidly  due  to  wavenumber  resolution:  due  to  the  coarse  grid  {d^  =  100)  Singer’s 
data  has  oj^  =  0.4;  the  LESS  data  has  a;+  =  0.8.  The  data  of  Farabee  and  Casarella^^ 
is  slightly  lower  than  the  data  of  Gravante  et  and  has  a  frequency  range  with  a 
constant  decay  rate  slightly  larger  than  The  reason  for  the  faster  decay  rate  is 
that  it  has  >  0.8.  Above  a;+  >  0.8  the  spectra  of  Farabee  and  Casarella^^  does  not 
decay  as  fast  as  the  simulations. 
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Table  2.  Wall  pressure  frequency  spectra  parameters.  Rco  =  Uo5/u,  where 
Up  is  the  channel  centerline  velocity  for  channel  flow  and  free 
stream  velocity  for  boundary  layer  flows. 


Source 

Rbq 

RCr 

d+ 

< 

LESS 

3200 

171 

44 

0.8 

DNSS 

3200 

180 

18 

2.1 

Singer 

25,500 

1000 

100 

0.4 

Farabee  and  Casarella^^ 

39,000 

1535 

44 

0.8 

Gravante  et  al.^^ 

40,000 

1504 

8 

5 
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Figure  4.  Comparison  of  (a)  outer-  and  (b)  inner-scaled  (5+  =  vjur)  point 

frequency  spectra;  —  :  LESS, - :  DNS5, - :  Singer. 

A  :  Farabee  and  Casarella^^,  and  □  :  Gravante  et 


kx  —  u)  spectra —  R1234 

In  this  subsection  we  show  the  two-dimensional  kx—oj  spectra  for  the  total  pressure 
from  all  regions  and  compare  it  to  that  from  the  literature  and  previous  simulations. 
Figure  5  shows  the  two-dimensional  spectra  from  LESS,  LES9  and  from  Choi  and  Moin.® 
(reduced  by  Bruce  Abraham  and  shown  in  Chang,  Abraham  and  PiomeUi^).  The  CM 
data  have  not  been  adjusted  for  the  change  in  normalization  and  only  the  wavenumber- 
frequency  range  corresponding  to  LESS  is  shown.  The  character  of  the  subconvective 
range  is  quite  different:  the  data  of  CM  has  much  higher  energy,  basically  flattening 
out,  as  has  been  noted  previously.  The  LESS  data  falls  off  rapidly  with  decreasing  |A:a;| 
with  the  energy  build-up  in  the  the  smallest  |A:a;|  modes;  the  LES9  data  is  flat  in  the 
smallest  |fca;|  modes. 

Figme  6  are  constant-o;  cuts  comparing  LESS,  LES9  and  CM.  The  CM  data  has 
a  convective  ridge  that  is  both  higher  and  wider  than  the  LES;  this  is  consistent  with 
previous  findings.  What  is  notable  is  that  in  the  subconvective  range  the  CM  data 
tends  toward  a  behavior  as  u  increases,  whereas  the  LESS  and  LES9  data  has  an 
increasingly  large  rolloff  in  this  range.  The  LES9  data  is  quite  different  than  the  LESS 
data  in  the  lowest  1^x1  modes:  it  is  flat,  not  having  the  energy  build-up  as  in  LESS. 

Partial  MS  and  TT  pressures  —  R1234 

Figme  7  shows  the  one-dimensional  spectra  of  total,  MS  and  TT  pressmes  due 
to  somces  over  the  entire  channel.  The  streamwise  and  spanwise  spectra  are  similar 
to  the  DNS  data  shown  in  Chang,  Piomelli  and  Blake^.  The  plots  show  that  the 
low  wavenumbers  and  frequencies  are  dominated  by  the  MS  pressme;  in  the  spanwise 
wavenumber  and  frequency  mid-range  both  TT  and  MS  contribute  significantly  to  the 
total  pressure;  at  high  spanwise  wavenumbers  the  TT  pressme  becomes  dominant  (see 
e.g.,  the  spanwise  spectra,  Figme  7(b)). 

Plots  of  the  two-dimensional  kx—oJ  spectra  of  the  total,  MS  and  TT  partial  pressmes 
due  to  somces  from  the  entire  channel  (R1234)  for  LESS  and  LES9  are  shown  in  Figmes 
8  and  9,  respectively.  Figme  8(a)  shows  that  the  total  pressme  has  an  energy  buildup  in 
the  lowest  wavenumbers.  Comparison  with  Figmes  8(b)  and  (c)  shows  that  this  buildup 
must  come  from  the  TT  pressme;  the  MS  pressure  falls  off  rapidly  in  the  lowest  [fca,! 
bins.  Both  the  TT  and  total  pressme  for  LES9,  Figme  9,  have  smooth  contoms  with 
no  buildup  in  the  lowest  |fcx|  bins.  The  MS  pressme  for  LES9  is  identical  to  that  shown 
for  LES8  in  Figme  8(b). 

Greater  insights  into  the  natme  of  the  subconvective  wavenumber  range  can  be 
obtained  from  constant  frequency  plots  of  streamwise  wavenumber  spectra.  Figure  10 
(LES8)  and  Figure  11  (LES9).  It  is  obvious  that  the  TT  pressme  is  the  somce  of 
the  subconvective  pressme  as  the  MS  pressure  goes  to  zero  as  kx  0.  At  the  lower 
frequencies,  Figme  10(a)  and  Figme  11(a),  the  subconvective  range  is  flat;  as  frequency 
increases  [Figmes  10(b)-(d)  (LES8)  and  Figmes  ll(b)-(d)  (LES9)],  the  energy  increases 
as  kx  0.  At  the  lower  frequencies  the  MS  and  TT  contributions  at  the  convective 
ridge  are  about  equal.  As  frequency  increases,  the  MS  becomes  the  dominant  pressme 
at  the  convective  ridge.  Though  the  wavenvunber  resolution  limits  om  ability  to  make  a 
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conclusive  statement,  it  appears  that  the  high  wavenumber  side  of  the  convective  ridge 
is  dominated  by  the  MS  pressure. 
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Figure  8.  Two-dimensional  contour  plots  of  A:*  —  a;  spectra  from  LESS  for 
all  regions  (R1234);  (a)  (b)  and  (c) 
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Figure  10.  Constant  frequency  cuts  of  fc®  —  u;  spectra  of  partial  pressures 
from  LESS  for  all  regions  (R1234);  (a)  u5/ur  =  52  (b)  u)5lur  = 
104  (c)  (jj5lur  =  156  and  (d)  ud/ur  =  204;  — —  :  total; - : 


Figure  11.  Constant  frequency  cuts  of  —  co  spectra  of  partial  presstures 
from  LES9  for  all  regions  (R1234);  (a)  uS/ur  —  44  (b)  uS/ur  = 


108  (c)  uSfur  =  156  and  (d)  (joSjur  =  204; 


total; 


Partial  pressures  from  regions  of  channel 

In  this  subsection  we  look  at  the  partial  pressures  from  four  regions  of  the  channel. 

One-dimensional  spectra  of  partial  pressures 

Figure  12  shows  the  one-dimensional  MS  pressure  spectra  for  all  regions  of  the 
channel  and  for  each  of  the  individual  regions.  The  streamwise  and  spanwise  spectra 
are  similar  to  the  DNS5  data  shown  in  Chang,  PiomeUi  and  Blake^.  The  frequency 
spectra  have  previously  not  been  shown.  The  figures  show  that  the  the  buffer  layer 
dominates  the  entire  wavenumber  and  frequency  ranges;  the  logarithmic  region  is  the 
next  most  dominant;  the  viscous  shear  layer  is  small  at  low  kx  —  u),  increasing  with 
kx  —  all  regions  except  the  viscous  shear  layer  are  important  at  low  kx  —  w. 

The  one-dimensional  TT  pressme  spectra  for  all  regions  of  the  channel  and  for 
each  of  the  individual  regions  from  LESS  are  shown  in  Figure  13.  The  plots  for  LES9 
are  virtually  identical  and  have  not  been  included.  It  shows  the  familiar  low  kx  —  ijO 
energy  buildup  shown  in  Chang,  Piomelli  and  BlakeL  The  low-frequency  energy  buildup 
indicates  that  the  partial  pressmes  have  high-energy  slowly-varying  components.  These 
have  been  observed  in  animations  of  the  TT  partial  pressmes. 

Two-dimensional  kx  —  oj  spectra  of  p£u*tial  pressures 

Figures  14  and  15  show  the  MS  spectra  for  sources  from  each  of  the  regions,  re¬ 
spectively.  It  is  apparent  that  with  increasing  distance  from  the  wall,  the  slope  of  the 
convective  ridge  in  kx  —  u)  space  {u)/{kx5)  increases;  the  width  of  the  convective  ridge 
decreases;  the  shape  of  the  convective  ridge  becomes  increasing  asymmetric;  the  energy 
in  the  higher  wavenumbers  and  frequencies  becomes  attenuated,  becoming  increasingly 
concentrated  in  low  kx  —  u).  The  TT  kx  —  oj  spectra  from  the  four  regions  are  shown 
in  Figures  16  and  17  for  LESS  and  Figures  18  and  19  for  LES9.  The  figures  show  that 
both  LESS  and  LES9  have  similar  energy  build-ups  in  the  low-wavenumbers.  However, 
as  shown  previously,  when  the  pressures  from  all  the  regions  are  added  together  LESS 
and  LES9  are  quite  different:  Figure  8(a)  and  (c)  shows  that  for  LESS  the  TT  partial 
pressure  (and  total,  too)  from  all  regions  retains  the  low-|A:a;|  energy  buildup.  On  the 
other  hand,  Figmre  9(a)  and  (b)  show  that  adding  the  contributions  from  all  regions 
results  in  the  low-|A:^|  energy  buildup  disappearing. 

Convection  velocities 

As  noted  above,  the  convection  velocities  computed  from  the  peak  of  the  convective 
ridge  can  characterize  the  kx  —  (jJ  behavior  of  the  peak.  These  were  computed  from 

v,{k,)  =  (22) 

and  are  shoAvn  in  Figure  20  for  the  MS  and  TT  partial  pressures  from  each  of  the 
regions.  Uc  increases  with  distance  from  the  wall,  from  about  lOUr  due  to  viscous 
shear  layer  sources  to  about  15ur  due  to  the  logarithmic  region  and  ISuj  due  to  the 


23 


upper  channel.  The  convection  velocities  from  the  MS  and  TT  partial  pressiures  for 
all  the  regions  shows  that  for  the  MS  partial  pressures,  the  low  wavenumbers  have  the 
same  value  as  the  logarithmic  region,  but  approach  the  buffer  layer  curve  at  the  higher 
wavenumbers.  The  peak  value  of  Uc  appears  to  be  15.3ur  at  kxS  =  3,  decreasing  to 
about  llttr  at  the  highest  wavenumbers. 
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Figure  12.  One-dimensional  spectra  of  tt^'®  from  LESS;  (a)  streamwise 
wavenumber  and  (b)  spanwise  wavenumber  and  (c)  frequency; 

- :  all  regions  (R1234),  - -  :  viscous  shear  layer(Rl), 

.  :  buffer  layer, - :  logarithmic  region,  and  —  : 

upper  channel. 


lOlogK(k) 


Figure  13.  One-dimensional  spectra  of  tt”’  from  LESS  (a)  streamwise 
wavenumber  and  (b)  spanwise  wavenumber  and  (c)  frequency. 

-  :  all  regions  (R1234),  -  :  viscous  shear  layer(Rl), 

.  :  buffer  layer(R2), - :  logarithmic  region(R3), 

and - :  upper  channel(R4). 
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Figure  15.  Two-dimensional  contour  plots  oi  kx  —  ijJ  spectra  of  froi 
LESS;  (a)  logarithmic  region  (R3)  and  (b)  upper  channel  (R4) 


Figure  17.  Two-dimensional  contour  plots  of  kx  —  u)  spectra  of  tt-'  ■'  froi 
LESS;  (a)  logarithmic  region  (R3)  and  (b)  upper  channel  (R4) 


Figure  19.  Two-dimensional  contoiu  plots  oi  kx  —  UJ  spectra  of  tt^  from 
LES9;  (a)  logarithmic  region  (R3)  and  (b)  upper  channel  (R4). 


Constant— a;  cuts  of  kx—u  spectra  of  partial  pressures 

Constant  frequency  cuts  of  the  kx  —  to  spectra  for  the  MS  pressures  from  all  the 
regions  compared  with  that  for  the  individual  regions  are  shown  in  Figure  21.  At  the 
lowest  frequency,  Figure  21(a),  the  low-^a;  side  of  the  convective  ridge  seems  to  be 
due  to  the  logarithmic  region;  the  peak  due  to  both  the  buffer  layer  and  logarithmic 
region;  the  high-fc^  side  due  to  the  buffer  layer.  As  frequency  increases,  Figure  21  (b- 
d),  the  importance  of  the  logarithmic  region  diminishes  and  the  buffer  layer  appears 
to  be  primarily  responsible  for  the  convective  ridge  energy.  In  the  two  lowest-o;  cases 
[Figure  21(a-b)]  the  subconvective  range  has  an  approximate  kl  behavior.  As  frequency 
increases,  the  subconvective  range  goes  to 

Constant  frequency  cuts  of  the  kx  —  oj  spectra  for  the  TT  pressures  from  LESS  from 
all  the  regions  compared  with  that  for  the  individual  regions  are  shown  in  Figure  22. 
Results  from  LES  were  almost  identical,  so  are  not  shown.  Similar  to  the  broad-band 
partial  pressures,  in  the  low-fc*  the  constant-o;  cuts  from  the  individual  regions  have 
higher  energy  than  that  due  to  all  the  regions  (R1234).  For  all  the  frequencies  the 
lowest-fcx  partial  pressures  from  the  buffer  layer  and  logarithmic  region  are  essentially 
equal,  diverging  as  kx  approaches  the  convective  ridge.  This  is  not  surprising,  in  light  of 
the  fact  that  the  one-dimensional  broad-band  wavenumber  spectra  for  the  buffer  layer 
and  logarithmic  region  are  very  close  in  the  the  lowest-fc^  (Figure  13).  With  regards  to 
the  partial  pressures  from  the  buffer  layer  and  logarithmic  region,  as  frequency  increases, 
the  relative  level  of  the  pressure  in  the  lowest  kx  decreases  with  respect  to  that  at  the 
convective  ridge;  this  points  out  that  the  \ow-kx  energy  build-up,  though  evident  at  all 
frequencies  is  mostly  concentrated  at  low-o;.  A  feature  of  the  partial  pressures  from  all 
the  regions  is  that  as  frequency  increases,  the  minimum  between  the  subconvective  and 
convective  ridge  peaks  becomes  more  pronormced  and  move  to  higher  kx- 

In  the  low-w  [Figure  22(a)],  the  \ow-kx  side  and  peak  of  the  convective  ridge  mimics 
the  behavior  of  the  logarithmic  region,  while  the  high-fcx  side  mimics  the  pressure  from 
the  buffer  layer.  Note  that  for  part  of  the  high-A:^  side  of  the  convective  ridge  the  partial 
pressure  from  the  buffer  layer  is  actually  higher  than  that  from  all  regions.  As  frequency 
increases,  the  convective  ridge  partial  pressure  due  to  the  logarithmic  region  decreases 
with  respect  to  that  from  the  buffer  layer;  thus  the  buffer  layer  partial  pressure  is  the 
dominant  source  of  convective  ridge  energy  there. 
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Comparison  to  semi-empirical  formulations 


The  incompressible  model  for  the  wall-pressure  spectra  of  Chase^®  is  given  by 


^c(K,a;)  = 


[Kl  +  (W)-2]"/2 


CtK^ 


Kl  +  jbS) 


-2 


-t-  (6(5)-2 


(23) 


where 

K2  -  ^ 

{hUr) 

There  are  four  constants  that  need  to  be  determined:  Ct  and  Cm,  which  determine  the 
amplitude  of  the  TT  and  MS  pressures,  respectively,  b  and  h.  The  convection  velocity, 
Uc  determines  at  which  value  of  kx  the  convective  ridge  appears.  It  has  been  assumed 
that  the  pressure  soiuce  layer  is  of  thickness  b5,  where  6  is  the  boundary  layer  thickness. 
(24)  takes  into  account  the  assumed  decorrelation  of  the  velocity-field  sources  as  they 
convect  downstream  at  Uc,  huj  is  a  characteristic  turbulence  velocity.  The  parameter 
h  controls,  to  some  extent,  the  width  of  the  convective  ridge.  Chase^®  recommended 
that  the  values  h  =  3,  Cxh  =  0.014,  =  0.466  and  b  —  0.75  be  used. 

Figure  23  compares  the  spectrum  of  Chase^®  to  the  total,  TT,  and  MS  spectra 
from  LES9  at  u5fur  =  156.  The  spectra  (23)  has  been  integrated  over  all  kz  to  be 
consistent  with  the  LES9  results.  The  results  of  Chase  are  shown  with  parameters 
which  lead  to  a  best  fit  of  the  subconvective  range  and  convective  ridge—  h  =  1, 
Crph  =  8  X  10”®,  CMh  =  5  X  10”^  and  6=1.  The  model  has  the  correct  wavenumber- 
white  [i.e.,4>{kx)  ~  k^]  behavior  in  the  subconvective  range,  and,  as  can  be  seen  in  (23) 
the  MS  component  has  a  k^  behavior  in  this  range  as  shown  e.^.,in  Figure  22.  The 
model  does  not  correctly  predict  the  width  of  the  convective  ridge —  this  could  be  due 
to  fact  that  it  considers  only  the  velocity-field  somces  in  the  logarithmic  region.  As 
shown  in  Figure  22  the  MS  pressure  from  the  logarithmic  region  is  indeed,  very  narrow- 
peaked,  and  the  full  width  of  the  convective  ridge  comes  from  sources  in  the  buffer 
layer.  The  rolloff  between  the  convective  ridge  and  subconvective  range  is  determined 
primarily  by  the  values  of  Cm  and  Ct,  and,  as  constants,  they  cannot  reflect  the  fact 
that  the  rolloff  varies  with  u  as  shown  in  Figure  22. 


CONCLUSIONS 

It  was  found  that  the  form  of  the  nonlinear  terms  makes  a  difference  in  the  TT 
pressure  for  the  entire  channel  in  the  subconvective  wavenumber  range.  With  the 
rotational  form  (LESS),  the  subconvective  pressures  have  high  energy  build-up  in  the 
lowest  wavenumbers,  whereas  with  the  skew-symmetric  form  (LES9),  there  is  no  such 
behavior.  When  the  subconvective  pressures  are  decomposed  into  those  from  the  four 
regions  of  the  channel,  the  effect  of  advective  scheme  is  not  apparent:  the  form  of  the 
nonlinear  terms  only  affects  the  pressure  over  the  entire  channel. 

It  was  shown  that  the  TT  pressure  is  exclusively  responsible  for  the  subconvective 
pressure  and  that  both  the  buffer  layer  and  logarithmic  region  are  the  main  contributors. 
The  MS  pressure  behaves  as  k^,  thus  going  to  zero  as  kx  0. 
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Figure  23.  Constant  frequency  cuts  of  —  a;  spectra  at  =  156  com¬ 
pared  with  Chase^^;  —  :  Chase; - :  total; - .  :  MS; 

A  :  TT. 


With  regards  to  the  pressures  at  the  convective  ridge,  the  one-dimensional  broad¬ 
band  spectra  shows  that  the  MS  pressure  is  dominant  in  the  low-  to  medium- wavenumbers 
and  frequencies,  but  in  the  higher  wavenumbers  the  TT  pressure  becomes  dominant. 
Constant  frequency  slices  of  the  two-dimensional  kx  —  oj  spectra  shows  that  in  the  lower 
frequencies  both  the  logarithmic  region  and  buffer  layer  are  responsible  for  the  MS  con¬ 
vective  pressures,  but  as  frequency  increase,  the  buffer  layer  becomes  more  dominant. 
The  same  is  true  for  the  TT  pressure. 

Spectra  predicted  with  the  semi-empirical  model  of  Chase^^  was  compared  with  LES 
results.  It  is  apparent  that  such  models  can  be  tuned  to  have  the  correct  behavior  in 
the  subconvective  range  and  height  of  the  convective  ridge.  However,  possibly  due  to 
lack  of  inclusion  of  buffer  layer  sources,  the  width  of  the  convective  ridge  is  underpre¬ 
dicted.  Also,  such  a  model  cannot  describe  the  frequency-dependent  rolloff  between  the 
convective  ridge  and  subconvective  wavenumber  range.  Parameterization  of  the  buffer 
layer  sources  is  necessary  for  a  fully-successful  wall-pressure  model. 
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